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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

TRANSONIC FLUTTER INVESTIGATION OF MODELS OF THE 
SWEPTBACK WING OF A FIGHTER AIRPLANE 


By Samuel L. Smith III and Robert W. Boswinkle, Jr. 


SUMMARY 


A transonic flutter investigation has been made of models of the 
wing of a current fighter airplane. The models were dynamically and 
elastically scaled in accordance with criteria which include a flutter 
safety margin. The wings had an aspect ratio of 3-42 and were swept 
back 4l.l° along the leading edge and 19-3° along the outer part of the 
trailing edge . A large trailing-edge fillet extended out to 50 percent 
of the semispan. The investigation was made in the Langley transonic 
blowdown tunnel and covered a Mach number range from 0-75 to I. 32 . 


The flutter boundary was located at simulated altitudes below sea 
level, the models being flutter free at altitudes above sea level. 
However, a region in which the models exhibited large responses to the 
turbulence of the tunnel stream extended to altitudes above sea level 
at supersonic Mach numbers . The significance with regard to the air- 
plane of the large responses of the models is not known. The flutter 
boundary shifted to higher altitudes but remained below sea level with 
the addition of 15 -percent-chord leading-edge extensions over the outer 
35 percent of the semispan. 


INTRODUCTION 


The flutter characteristics of the wing of a current fighter air- 
plane have been under study. The wing is swept back 41.1° along the 
leading edge and 19*3° along the outer part of the trailing edge. A 
large trailing-edge fillet extends out to 50 percent of the semispan. 
Calculations indicated that flutter would result at transonic speeds 
at sea level if the stiffness were reduced only slightly. Experimental 
data on similar wings (refs . 1 to 4) indicated that possibly a sufficient 


stiffness margin existed; however, it was felt that the wing in question 
was sufficiently different from those of the references to warrant a 
separate experimental study. 

The investigation was made in the Langley transonic blowdown tun- 
nel with models which were dynamically and elastically scaled in accord- 
ance with criteria which include a flutter safety margin. The wing 
spar was cantilever-mounted inboard of the wing root and the tests were 
made at Mach numbers from 0.75 to 1.32 and at simulated altitudes 
extending to below sea level. The effect of installing a 15-percent- 
chord leading-edge extension over the outer 35 percent of the semispan 
was also investigated. 


SYMBOLS 


b typical wing semichord, ft 

c local streamwise chord, ft 


l 


m 


length scale factor, 


Typical model length 

Corresponding airplane length 


mass scale factor. 


Typical model mass 

Corresponding airplane mass 


m' mass of exposed panel, slugs 


M Mach number 

q dynamic pressure, lb/sq ft 

s value of y at wing tip 

t time scale factor. 

Time required for tunnel airstream to move 1 model chord length 
Time required for airplane to move 1 airplane chord length 

T static temperature , °R 


v 


t fo ° 2 ^ 


V 


velocity, ft/sec 
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V reduced velocity based on a representative natural frequency, 

V/bo^ 

y distance from wing root measured perpendicular to wing root, ft 

X,Y streamwise and spanwise coordinates, respectively, defined 

in figure 4 

T) stiffness reduction factor used to provide margin of safety 

in application of model flutter-test results to the airplane 

p mass ratio., m'/pv 

p static air density, slugs/cu ft 

o>j_ representative natural frequency, radians/sec 

Subscripts : 

A airplane 

M model 


MODELS 


Geometry 

The models were 5.125-percent-size versions of the wings of a 
current fighter airplane. The wing models had an aspect ratio of 3.42 
and were swept back 4l.l° along the leading edge and 19-3° along the 
outer part of the trailing edge. A large fillet at the trailing edge 
extended out to 50 percent of the semispan. A sketch of the model is 
given in figure 1 and some of the more important geometric properties 
are listed in table I. The fact that the plan-form aspect ratio is 
twice the exposed-panel aspect ratio (table I) is coincidental. 

Because of damage to the models at flutter, six models were 
required in the investigation. Three models (designated wings 1 to 3) 
were without leading-edge chord-extensions and were intended to be 
identical. The other three models (designated wings 4 to 6) had leading- 
edge chord-extensions and were intended to be identical. In addition, 
the only intended differences between the two sets of models were 
differences caused by the addition of the leading-edge chord-extensions. 
Small differences between models 1 to 3 and also between models 4 to 6 
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did exist, as evidenced by the measured natural vibration frequencies 
and node lines (presented in the section entitled "Physical Properties"). 

The chord-extensions were over the outer 35 percent of the semi- 
span and increased the local wing chords by 15 percent . A model with 
leading-edge chord-extensions is shown mounted in the fuselage mounting 
block in figure 2. (As shown in figure 2, the wings were painted at 
intervals along the leading edge to aid in observing the motion of the 
models during the flutter runs . ) The wings without leading-edge chord- 
extensions had a small amount of positive camber and the leading-edge 
chord-extensions of models k to 6 accentuated the camber. 


Scaling 

The nondimens ional mass and stiffness distributions were required 
to be the same for the model as for the airplane. The mass and stiff- 
ness levels for the model were obtained by specifying the scale factors 
for the fundamental quantities involved: length, mass, and time. 

The size of the model was limited by tunnel-wall-interference 
effects, and on the basis of past experience the length scale factor 
was chosen to be 


l = 0.03125 


( 1 ) 


The mass scale factor was obtained from a requirement that the 
mass ratio p should be the same for the model as for the airplane, 
which results in 


_ P M 7 3 

m = — l 

P A 


( 2 ) 


In order to locate the simulated sea level near the middle of the tunnel 
density range available at a Mach number of 1, the density ratio was 
chosen to be PyijPp^ = 1*97* This location of simulated sea level allows 

altitudes below sea level to be obtained and flutter margins to be 
indicated for cases where flutter does not occur above sea level. 

The time scale factor was obtained from a requirement that the 
reduced velocity V should be the same for the model as for the air- 
plane, which results in 
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t = 



l 


Since the Mach number is the same for the model as for the airplane, 
the time scale factor may be written 



(3) 


The static temperature for the airplane T A is a function of altitude 

only, and for sea level it was taken to be 519° R. However, in the 
tunnel during a run, the temperature continually drops as air is expended 
from the reservoir and the temperatures obtained at the various flutter 
points during an investigation are different. A study of previous flut- 
ter data indicated that 4o8° R was near the average value of the static 
temperature that would be expected during the present runs, and this 
value was used to obtain the temperature ratio used in the scaling: 

t m/ t a = 0 ‘ 786 * 

A list of the pertinent wing and flow quantities and the design 
scale factors used is given in table II. It may be noted that the 
factor q is used in the scale factors for some of the quantities 
listed. The factor q has the value 0.7 6 and occurs because the 
stiffnesses of the models were made percent of those which would 
result from application of the scale factors as specified (eqs. (l), 

(2), and (3))* The purpose of reducing the model stiffnesses was to 
provide a margin of safety in the application of the model flutter- 
test results to the airplane. Thus the design reduced velocity for the 
model is equal, not to that of the airplane, but to that of an airplane 
having stiffnesses 76 percent of those of the actual airplane. 

The dynamic pressure and Mach number are quantities which are 
controllable during a run, whereas the temperature is not. If the 
dynamic pressure and Mach number are considered to be fixed and a 
static temperature different from the design value is obtained, both 
the density and velocity will be different from the values considered 
in the scaling. The density and velocity changes result, respectively, 
in values of mass ratio and reduced velocity different from the design 
values . However, a combination of reduced velocity and mass ratio 
which can be expressed in terms of the dynamic pressure 
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is independent of the temperature, and this combination is exactly simu- 
lated in the runs by the expedient of interpreting the simulated alti- 
tude in terms of dynamic pressure. Thus, the scale factor in table II 
for dynamic pressure is used to convert the dynamic pressure for the 
airplane at any altitude and Mach number to the dynamic pressure for 
the model at the same altitude and Mach number. The dynamic pressure 
for the airplane is assumed to be that calculated by use of the ICAO 
standard atmosphere (ref. 5 )* It may te noted that, for a given alti- 
tude, q/M^ is a constant. 

The effect of not having the mass ratio and reduced velocity of 
the models exactly equal to those of the airplane is believed to be 
negligible in the present investigation. Experience with a wide variety 
of flutter models has indicated that, at least within the operational 
limits of the tunnel, flutter at a given Mach number tends to occur at 
a constant value of dynamic pressure regardless of the individual values 
of density and velocity. 


Construction 

The construction of the models is indicated in figure 1. The main 
spar was made of aluminum alloy, and aluminum-alloy ribs having U-shaped 
cross sections were welded to the main spar. The leading and trailing 
edges were of pine. Balsa was used to fill the wing to contour. Lead 
weights were placed in the wing at various locations and the wings were 
wrapped with silk cloth and painted. Each wing panel was instrumented 
with strain gages on the main spar near the root . The main spar was 
clamped inboard of the root, as shown in figure 1 , and thus allowed 
some root flexibility. The mounting block shown in figure 2 was made 
of aluminum alloy. 


Physical Properties 

The first several natural cantilever frequencies and node lines of 
each of the six wings are given in figure 3 . In obtaining the data an 
electromagnetic shaker was used to excite each panel separately . The 
shaker stem acted on the extended wing spars at the locations indicated 
by x in figure 3 and- the spars were clamped as indicated in figure 1 . 
The positions of the node lines were indicated by salt crystals sprinkled 
on the wings . 
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The right panel of model 2, which survived the flutter tests 
undamaged, was used to obtain the flexibility influence coefficients. 
Influence coefficients were obtained at 22 stations (fig. *0 on the wing 
by the method described in reference 6. The influence-coefficient matrix 
is given in table III. This matrix has been made symmetrical in table TV 
by taking the average of each pair of coefficients symmetric to the diag- 
onal. The deviation of the coefficients in table III from the average 
values in table IV gives some indication of the accuracy of the measure- 
ments. Only 2.6 percent of the coefficients deviate more than 2 percent, 
and the greatest deviation is 3 .6 percent. 

The right panel of model 2 was cut into strips and the center of 
gravity, mass, and moment of inertia about the center of gravity of 
each strip were measured. The data are given in figure 5. Each strip 
was then cut as shown in figure 4 so that each section corresponded to 
one of the influence coefficient stations . The mass and center of 
gravity of each section were measured and the values are listed in fig- 
ure 4 . The masses given in figures 4 and 5 Tor the sections and strips 
include an allowance for the material lost in the saw cuts . 


APPARATUS AND TESTS 


The investigation was made in the Langley transonic blowdown 
tunnel, which has a slotted test section. The test section is octagonal 
in cross section and measures 26 inches between flats. During the oper- 
ation of the tunnel, a preselected Mach number is set by means of a 
variable orifice downstream of the test section, and this Mach number 
is held approximately constant (after the orifice is choked) while the 
stagnation pressure, and thus the density, is increased. The static 
density range is approximately 0.001 to 0.012 slug per cubic foot, and 
Mach numbers from subsonic values to a maximum of about 1.4 may be 
obtained. Because of the expansion of the air in the reservoir during 
a run, the stagnation temperature continually decreases, and therefore 
the test-section velocity is not uniquely defined by the Mach number. 
Additional details of the tunnel are contained in reference 1. Excel- 
lent agreement between flutter data obtained in the tunnel and in free 
air has been observed (ref. 7)* 

In the investigation, each model was cantilever-mounted in the 
mounting block shown in figure 2. The mounting block was fitted into 
a sting in such a way as to form a fuselage 3 inches in diameter which 
extended upstream into the subsonic flow region of the tunnel. This 
arrangement prevented the formation of shock waves from the fuselage 
nose which might reflect from the tunnel walls onto the model. A 
sketch of the model mounted on the sting and installed in the tunnel 
is shown in figure 6. The sting and model weighed approximately 
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290 pounds and the system had a fundamental bending frequency of about 
15 cycles per second. 

Wire strain gages were mounted on the wing spars near the root and 
were oriented so as to indicate model deflections about two different 
axes. The strain-gage signals, the tunnel stagnation and static pres- 
sures, and the stagnation temperature were recorded on a recording 
oscillograph. The strain-gage signals were used to indicate the start 
of flutter and the flutter frequency. High-speed motion pictures were 
made during some of the runs. 

The wings without leading-edge chord-extensions were tested at 
zero angle of attack. The wings with leading-edge chord-extensions 
were tested at - 2 ° angle of attack in an attempt to reduce the static 
loads . 


RESULTS AND DISCUSSION 


Presentation of Data 

The results of the investigation are given in table V(a) for the 
wings without leading-edge chord-extensions and in table V(b) for the 
wings with leading-edge chord-extensions. The dynamic pressure at the 
various test points is plotted as a function of Mach number in figure 7 
for the wings without leading-edge chord-extensions and in figure 8 
for the wings with leading-edge chord-extensions. Lines of constant 
simulated altitude are also indicated in figures 7 and 8 . 

Each circle symbol in figures 7 and 8 indicates the point of the 
start of definite flutter and each square symbol indicates the point 
of the maximum dynamic pressure attained during a run without obtaining 
flutter. A dashed line below a symbol defines a low-damping condition. 
In the low-damping condition, the strain-gage records and the motion 
pictures indicated periods of nearly sinusoidal, lowly damped oscilla- 
tions. The point for the beginning of low damping in each run was 
indefinite and was somewhat arbitrarily chosen. On the other hand, the 
point for the beginning of flutter in each run in which flutter was 
obtained was definite and was characterized by rapidly diverging oscil- 
lations . The low-damping region is indicated for the wings without 
leading-edge chord-extensions in figure 7 8 y dotted shading. 

The response frequencies of the wings are indicated near most of 
the data points in figures 7 and 8 . The response frequency for no- 
flutter or low-damping points was taken as the predominant oscillation 
frequency of the models; at flutter, of course, the flutter frequency 
is listed. 
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The flutter mode for both configurations investigated involved 
bending and torsion of the wing with some rotation in pitch at the wing 
root. The rotation in pitch of the wing root was possible because, 
as previously noted, the main spar was clamped inboard of the root. 

A typical oscillograph record showing the strain-gage traces during 
low damping and flutter is given in figure 


Interpretation of Results 

As stated in the section entitled "Scaling," the stiffnesses of 
models were 76 percent of the scaled airplane stiffnesses. The simulated 
altitudes indicated in figures 7 and 8 are thus to be interpreted as 
altitudes which, if cleared by the model, could be reached with a 
52-percent (1/0.76 = 1-32) margin of safety in stiffness by the airplane. 
This statement assumes, of course, that in all other respects the model 
exactly simulates the airplane. 

An alternate interpretation of the results arises from the fact 
that for most configurations the dynamic pressure required for flutter 
varies, to a first approximation, directly with the stiffness level. 

Thus, a flutter point obtained with the model indicates that the air- 
plane will flutter at the same Mach number at a simulated altitude 
corresponding to a dynamic pressure 32 percent higher than that for the 
model . 


Wings Without Leading-Edge Chord-Extensions 

The transonic flutter boundary for the models of the wing without 
leading-edge chord-extensions is located at altitudes below sea level 
(fig. 7) • The dynamic pressure for flutter is indicated to be a mini- 
mum at a Mach number of about 0.87- The low-damping region extends at 
supersonic Mach numbers to altitudes above sea level. With regard to 
the airplane, the significance of the low damping obtained with the 
models is not known. Photographs of the wings without leading-edge 
chord-extensions after flutter are given in figures 10(a) to 10(c). 


Wings With Leading-Edge Chord-Extensions 

Because of various data-recording difficulties, the flutter points 
at the three lowest Mach numbers for the wings with leading-edge exten- 
sions (fig. 8) are known only to an estimated accuracy of ±100 lb/sq ft 
for dynamic pressure and ±0.03 for Mach number. However, the shape of 
the transonic flutter boundary is shown to be similar to that for the 
wings without leading-edge chord-extensions (fig. 7)* Although the 
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flutter boundary shifted to higher altitudes with the addition of the 
leading-edge chord-extensions , no flutter was obtained at altitudes 
above sea level. 

Low damping preceded the flutter points at the lowest Mach numbers, 
but the location of these points could not be ascertained and they are 
omitted in figure 8 and table V(b) . A photograph of one of the wings 
with leading-edge chord-extensions after flutter is given in figure 10(d) . 


CONCLUSIONS 


The transonic flutter characteristics of models of the sweptback 
wing of a current fighter airplane have been studied in the Langley 
transonic blowdown tunnel. The models were dynamically and elastically 
scaled in accordance with criteria which include a flutter safety margin. 
The scaling was such that if at a given Mach number a certain altitude 
is cleared by the model, that Mach number and altitude could be reached 
with a 52 percent margin of safety in stiffness by the airplane. The 
following results were obtained: 

1. Although the flutter boundary for the wings without leading- 
edge chord-extensions was located at altitudes below sea level, a region 
of lowly damped oscillations that extended to altitudes above sea level 
was obtained at supersonic Mach numbers. 

2. With the addition of 15-percent-chord leading-edge extensions 
over the outer 55 percent of the semispan, the flutter boundary shifted 
to higher altitudes but remained below sea level. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., December 20, 1957 • 
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TABLE I.- GEOMETRY OF MODELS WITHOUT 
LEADING-EDGE CHORD-EXTENSIONS 


Streamwise airfoil section, tip Modified NACA 65 AOO 6 

Streamwise airfoil section, root Modified NACA 65 AOO 7 

Leading-edge sweepback, deg 4l.l 

Trailing-edge sweepback, deg 19*3 

Span, ft 1.252 

Plan-form area based on extension of panels 

to model center line, sq ft 0.4582 

Plan-form aspect ratio based on extension of 

panels to model center line 3*42 

Fuselage diameter, ft 0.250 

Exposed -panel span, ft 0.498 

Exposed -panel area, sq ft 0.1453 

Exposed-panel aspect ratio 1*71 




TABLE II.- DESIGN SCALE FACTORS OF PERTINENT 
WING AND FLOW QUANTITIES 


^ = 1.97; ^ = 0.786; t ) = 0.76 

P A t A 


Quantity 

Design scale factor 

Symbolical 

Numerical 

Fundainental quantities : 



Length 

1 

1 

0.03125 

Mass 

_ P M ,5 
m = l 

0.6012 X 10"^ 

Time 

P A 

/T \"l/ 2 

t = ( M 1 

\\) 

0.03525 

Derived quantities : 



Stream velocity 

It" 1 

0.886 

Stream dynamic pressure 

B 

C-J 

1 

cT 

1 

ro 

1.548 

Moment of inertia 

ml 2 

O.587 x io-7 

Flexibility influence coefficients . 

TfW 

27.195 

Natural vibration frequencies . . . 

Wvi 

24.73 

Bending and torsional stiffnesses . 

rjl^mt” 2 

1'. 122 X 10" 6 
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Figure 2.- Photograph of wing with leading-edge chord-extensions in mounting block. (Wings 

painted at intervals along the leading edge.) 





3.- Measured natural vibration frequencies and node lines. Symbol x indicates shaker 
location. Numbers beside node lines indicate frequencies in cycles per second. 




Figure 4 -.- Sketch of wing without leading -edge chord -extensions, showing influence -coefficient 

station and center of gravity of various segments. 
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O Definite start of flutter 
□ Ifexiroun dynamic pressure, no flutter 

Low-damping condition 

Low-damping region 
Numbers beside data points indicate 
response frequencies in cpe 

Simulated altitude, ft = -10,00s7 ~T 
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Figure 7*- Transonic flutter characteristics of wings without leading 

edge chord-extensions. 
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O Definite start of flutter 

'O Definite start of flutter; location 

of flutter point estimated 
Loi*-damping condition 
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Figure 8.- Comparison of transonic flutter characteristics of wings with 
and without leading -edge chord-extensions. (For runs It, 15, and 17 
the accuracy of the data is less than that for the other runs, and 
although low-damping conditions preceded flutter, they are not 
indicated here.) 






Figure 9 .- A typical oscillograph record (run 4, wing l) 
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(c) Wing 3, front view. (d) Wing 5, top view. 

L-57-5501 

Figure 10.- Photographs of damaged models. (Wings were painted at intervals along 

the leading edge.) 






